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ABSTRACT 

The galactic black hole candidate Cygnus X-1 is observed to be in one of two X-ray 
spectral states; either the low/hard (low soft X-ray flux and a flat power law tail) or 
high/soft (high blackbody-like soft X-ray flux and a steep power law tail) state. The 
physical origin of these two states is unclear. We present here a model of an ionized 
accretion disc, the spectrum of which is blurred by relativistic effects, and fit it to 
the ASCA, Ginga and EXOSAT data of Cygnus X-1 in both spectral states. We 
confirm that relativistic blurring provides a much better fit to the low/hard state data 
and, contrary to some previous results, find the data of both states to be consistent 
with an ionized thin accretion disc with a reflected fraction of unity extending to the 
innermost stable circular orbit around the black hole. Our model is an alternative to 
those which, in the low/hard state, require the accretion disc to be truncated at a few 
tens of Schwarzschild radii, within which there is a Thomson-thin, hot accretion flow. 
We suggest a mechanism that may cause the changes in spectral state. 

Key words: accretion, accretion discs — binaries: close — black hole physics — 
stars: individual (Cygnus X-1) — X-ray: general — X-ray: stars 



> 1 INTRODUCTION 



Cygnus X-1 (Cyg X-1) is one of the brightest and best stud- 
ied X-ray sources in the sky. It was the first source to be 
identified with a binary system in which the X-ray emis- 
sion arises through accretion onto a compact object. (See 
e.g. Tanaka & Lewin (1995) for a review of black hole bina- 
ries). It is at a distance of 2 kpc (Massey et al. 1995) and 
consists of a supergiant star and a compact object, with an 
orbital period of 5.6 days. The mass of the unseen compan- 
ion is significantly greater than 5Mq (Dolan 1992) suggest- 
ing that it is a black hole. The high energy emission results 
from focused wind accretion from the supergiant (Gies & 
Bolton 1986) that is almost filling its Roche lobe. The X- 
ray fiux varies on all timescales from milliseconds to months 
although most of its time is spent in the so-called low or 
hard state (the 'low/hard' state) characterized by a low soft 
X-ray fiux and strong hard X-ray fiux. Occasionally Cyg X-1 
makes a transition from the low/hard state to the high/soft 
state which is characterized by a larger soft X-ray fiux (by 
a factor of ~ 10) and weaker hard X-ray flux. 

The spectrum in the low/hard state above a few keV 
is typically a power law of photon index T ~ 1.6-1.8 and a 
Compton reflection component above about 10 keV (Done 



et al. 1992; Ebisawa et al 1996) with a high-energy cut-off 
above ~ 100 keV. There is evidence of complicated iron ab- 
sorption and emission features around 6-7 keV (Barr, White 
& Page 1985; MarshaU et al. 1993; Ebisawa et al. 1996). 

In the high/soft state the spectrum is dominated by a 
soft blackbody-like component peaking around 1 keV with 
a power law tail of photon index T ~ 2.5 (Gierlinski et al. 
1999). 

It has been suggested that the different spectral states 
are due to a qualitative change in the accretion fiow (see e.g. 
Gierlinski et al. 1997; Esin et al. 1998; Poutanen & Coppi 
1998; Done & Zycki 1999). In the high/soft state there is a 
'standard' accretion disc extending in to the innermost sta- 
ble circular orbit around the black hole, and in the low/hard 
state the region within a few 10s of Schwarzschild radii be- 
comes a geometrically thick, optically thin, radiatively in- 
efficient advective flow (see also Nowak et al. 1999). The 
evidence for this is that in the low/hard state the reflected 
fraction appears to be low, which is consistent with the in- 
nermost parts of the accretion disc being missing, or at least 
Thomson thin. The change between spectral states occurs 
with a small change in bolometric luminosity (~ 50 per cent) 
and no change in 1.3 — 200 keV luminosity to within ~ 15 
per cent (Zhang et al. 1997). 
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Nov 
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Low /hard 
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Nov 
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1991 


Jun 6 


Ginga no. 2 


Low /hard 


b 


1996 


May 


SIS 


High/soft 


EOS 


1984 


Jul 9 


EXOSAT GSPC 


Low /hard 


E09 


1985 


Sep 14 


EXOSAT GSPC 


Low/ hard 


ElO 


1985 


Aug 12 


£;xo5'yir GSPC 


Low/ hard 


E13 


1984 


Nov 2 


.EXO^yir GSPC 


Low/ hard 


E14 


1984 


Nov 2 


EXOSAT GSPC 


Low /hard 



Table 1. Data sets used. The data set numbers 1 and 3-10 cor- 
respond to those of Ebisawa et al. (1996), and those prefixed by 
'E' correspond to those of Done et al. (1992). 



In this paper we investigate an alternative model in 
which the low inferred reflected fraction is a result of the 
ionization state of the accretion disc. A highly ionized disc 
may be an almost perfect reflector with weak reflection fea- 
tures (as compared to cold reflection) and hence may lead 
to an underestimate of the reflected fraction (Ross, Fabian 
& Young 1999; Ross, Fabian & Brandt 1996; Ross & Fabian 
1993; Nayakshin, Kazanas & Kallman 2000). The attraction 
of this model is that it does not require the nature of the 
accretion flow to change qualitatively between the spectral 
states; a change which is observed to occur with little varia- 
tion in bolometric luminosity as observed. It is just a change 
in density at the surface of the disc. 

Previous conclusions that the disc cannot be signifi- 
cantly ionized have been based on the absence of a strong 
photoelectric absorption edge at 8 keV. In the case of an ion- 
ized disc absorption and emission may occur at significant 
Thomson depths (tt ^ 1) and sharp features are blurred 
due to Compton scattering. See e.g. Fig. 7 of Ross, Fabian 
& Young (1999). 

In view of the extremely rapid continuum variability 
of Cyg X-1 we may also expect its spectrum to change 
rapidly. The spectra that we fit are necessarily time averaged 
since a typical observation will last for millions of dynamical 
timescales. 



2 DATA AND MODEL 
2.1 Data 

We use ASCA SIS, ASCA GIS, EXOSAT and Gmga data 
of Cyg X-1 in the low/hard state and ASCA SIS data of 
Cyg X-1 in the high/soft state. We do not use the ASCA 
93 Oct data of the low/hard state since this was taken when 
the satellite was in daytime (Ebisawa et al. 1996). The data 
used are summarized in Table ll 




E (keV) 

Figure 1. Reflection spectra from a slab of gas illuminated by an 
X-ray power law of photon index F = 1.6, solar iron abundance 
and various values of the ionization parameter ^. The lowest solid 
line is for log^ = 1, and log^ increases upwards to log^ = 5 in 
steps of 0.5. The dotted line is a representative illuminating power 
law, in this instance the log ^ = 3 case. 



2.2 Rest frame reflection spectra 

The accretion disc model is based on the X-ray refiection 
spectra from ionized slabs of gas calculated by Ross, Fabian 
& Young (1999) and Ross & Fabian (1993). The reflection 
spectrum resulting from the illumination of a slab of gas by 
an X-ray power law of photon index F depends primarily on 
the ionization parameter ^, defined as the ratio of flux and 
density, 

^ _ 47rF. 

where is the illuminating X-ray flux (0.01-100 keV) and 
wh is the hydrogen number density. Fig. ^ shows reflection 
spectra calculated for solar iron abundance, a photon index 
of 1.6 and a range of ionization parameters, log^ = 1 — 5 in 
steps of 0.5. 

For very low ionization parameters the most prominent 
feature of the reflected spectrum is the iron fluorescence line 
at 6.4 keV. In neutral gas both oxygen and iron are strong 
absorbers but the absorption above and below the iron edge 
is so strong that the change in observed flux across the edge 
is small, particularly if the primary continuum is also ob- 
served. This absorption also tends to harden the spectrum 
as compared to the illuminating power-law. 

As the X-ray illumination of the slab becomes more in- 
tense the matter becomes more highly ionized. This reduces 
the effects of absorption progressively from lower to higher 
energies. An important signature of ionized reflection is a 
strong iron edge. As much of the oxygen becomes completely 
ionized the edge becomes more prominent as the absorption 
below it is reduced. As the outer layers of the slab become 
completely ionized fluorescent photons are produced deeper 
in the slab and are Compton scattered on leaving the slab 
resulting in broad and blended emission and absorption fea- 
tures, e.g. the case of ^ = 10^ in Fig. ^. The narrow line cores 
are from those line photons that escape without scattering or 
are generated closer to the surface. For an extremely highly 
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ionization parameter ^ = 10^ the reflection spectrum is al- 
most indistinguishable from the illuminating continuum. 

2.3 Simple accretion disc model 

Initially we consider a model of a non-relativistic accretion 
disc with a single ionization parameter 5. This consists of two 
components; the primary X-ray source which we assume to 
be a power law of photon index F and the reflected spec- 
trum corresponding to some fraction of these X-rays being 
Compton scattered by the disc into our line of sight. 

We have computed a grid of reflection spectra for a 
range of photon indices, ionization parameters and iron 
abundances, and subsequently converted to an xspec 'table- 
model'. XSPEC is the spectral fitting code commonly used by 
X-ray astronomers. The reflection spectra in the tablemodel 
have been scaled so that, in xspec, a reflection spectrum 
plus power law of the same photon index has a correspond- 
ing reflected fraction equal to the ratio of the normalization 
of the reflected spectrum to the power law. We have consid- 
ered ionization parameters in the range ^ = 10^-10^ (up to 
^ = 10® in the high/soft case) and iron abundances in the 
range 1-3 x solar. 

2.4 Relativistic blurring 

The accretion disc is deep in the potential well of the black 
hole and subject to strong Doppler, transverse Doppler and 
gravitational redshifting effects. These will smear out the 
sharp features in the spectrum and blend the emission and 
absorption features together. The detection of such an ef- 
fect in the spectrum of Cyg X-1 in the low/hard state has 
been reported by Done & Zycki (1999) but with quite dif- 
ferent parameters to those reported here. We compute this 
in XSPEC by convolving the spectrum with a blurring kernel 
taken from the diskline model for an accretion disc around 
a Schwarzschild black hole (Fabian et al. 1989) . Special care 
needs to be taken when using this approach in xspec since it 
requires the model to be evaluated outside the energy range 
of the data. This is discussed in Appendix A. 

The parameters of the blurring model define the inner 
and outer radii of our accretion disc. The outer radius is fixed 
at lOOOrg = lOOOGM/c^, and the inner radius rin must be > 
6rg [Qrg — (SGM/c^ = 37?s, where Rs is the Schwarzschild 
radius, and is the location of the innermost stable circular 
orbit around a Schqarzschild black hole). The precise value 
of the outer radius heis little impact on the results we present 
here. The 'radial emissivity profile' para meter required for 
the diskline kernel is set to (1 - \/&/rg)/rl. 

2.5 A more realistic disc model 

We shall find that some of the data sets require the use of 
a more sophisticated model. In the future it will be useful 
to construct a proper model that accounts for the radial 
variation of ionization parameter in more detail, but for the 
time being we use the prescription described below. 

The accretion disc is divided into two zones, an 'inner' 
region which will have a relatively high ionization parame- 
ter, and an 'outer' region with a lower ionization parameter. 
The inner and outer regions are convolved with diskline 



kernels extending between 6r-g — > rtrans and rtrans — » lOOOrg 
respectively, where rtrans represents the radius of the tran- 
sition between high and low ionization zones. The relative 
contribution to the reflection from regions r < rtrans and 
r > rtrans are calculated to agree with the Page & Thorne 
(1974) radial emissivity law. From a practical point of view, 
fltting the data using xspec, this achieved iteratively, i) flt- 
ting the data, ii) fixing rtrans at the appropriate value given 
the relative contribution to the refiection from the two re- 
gions, iii) repeating until an self-consistent fit is achieved. 
Very few itierations were required. The error in determining 
rtrans is qulte small as an uncertainty of ~ 10 per cent in 
the ratio of fiuxes would correspond to an error of ±2rg in 
rtrans- This reprcscuts a realistic model of an accretion disc 
with a sharp transition from high to low ionization state. 
Whether such an abrupt transition is required will be inves- 
tigated later. 



3 FITS TO ASCA SIS DATA 

In this section we consider various fits to the ASCA SIS data 
of Cyg X-1 in both the high/soft and low/hard states. We 
are mainly interested in fitting the data between 4-10 keV 
since this is where most of the ionized refiection features 
are expected to lie. It is also uncertain as to whether the 
primary X-ray continuum is actually a power-law over the 
entire ASCA energy range. 



3.1 Simple disc model 

3.1.1 Low /hard state 

We begin by fitting a non-relativistically blurred ionized disc 
model to the 4-10 keV ASCA SIS data of the low/hard state. 
In these fits the iron abundance has been fixed at twice solar, 
the best fitting value suggested by Done & Zycki (1999). 
Such a model is not able to give a good fit to the data but 
it is interesting to consider a plot of against ^ shown 
in Fig. ^. The solid curve shows the quality of fit with the 
refiected fraction / frozen at 1, i.e. with the disc subtending 
2n sr of the sky as seen by the X-ray source. There are 
two local minima, one at ^ = 10'^ and the other around 
^ = 10*. The best fitting solution obtained will depend upon 
the initial value of ^, the separatrix being at approximately 
^ = 10^ . The horizontal dotted line shows the value for a 
refiected fraction of 0, and the dashed line shows the values 
with the reflected fraction free parameter. 

In the case of a low ionization parameter ^ = 10 a 
flt with a reflected fraction of 1 is signiflcantly worse than 
that of just a power law. A much smaller reflected fraction 
/ = 0.07 provides a better fit. For an ionization parameter 
of ^ = 10^ at the first minima in we again find a low 
refiected fraction is preferred, with / = 0.55. At the second 
minima ^ = 10^'*, however, we find the best fitting reflected 
fraction to be / = 1.07. Table | summarizes the results of 
these provisional fits. The values of the photon index F are 
consistent with the results of Gierliriski et al. (1997). 

To summarize, we have found two distinct solutions, 
one preferring a low ionization parameter and low refiected 
fraction, the other preferring high ionization parameter and 



© 0000 RAS, MNRAS 000, 000-000 



4 Young et al. 




log f 

Figure 2. (Low/hard). Results of fitting the ionized disc model, 
without relativistic blurring or additional spectral components, 
to the 93 Nov ASCA data of Cyg X-1 in the low/hard state. The 
solid curve shows the reduced value for a refiected fraction 
of 1, the dotted line for a refiected fraction of and the dashed 
line for the refiected fraction being a free parameter. When the 
reflected fraction is a free parameter it has been fitted for at each 
value of ^. 



Data set 


State 


log? 


r 




xVdof 


5 


low/hard 


1.0" 


1.57 


0.07 


128/76 


5 


low/hard 


2.0" 


1.65 


0.55 


78/76 


5 


low/hard 


3.8" 


1.68 


1.07 


75/76 


5 


low/hard 


5.0" 


1.54 


0.00 


178/76 


b 


high / soft 


1.0" 


2.36 


0.17 


57/37 


b 


high / soft 


2.4" 


2.47 


1.14 


33/37 


b 


high/soft 


5.1" 


2.32 


1.13 


37/37 


b 


high / soft 


6.0" 


2.21 


162 


70/37 



" frozen 

^ refiected fraction (2it sr) 

Table 2. Fit statistics for an ionized accretion disc model, with- 
out relativistic blurring or additional spectral components. Al- 
though some of these fits appear to be statistically acceptable 
some show systematic residuals at the expected energies of line 
and edge features. 

a high reflected fraction. We shaU investigate both of these 
with a more reahstic physical model. 

3.1.2 High/soft state 

Similar results are obtained for the high/soft state, where 
we have considered ionization parameters up to ^ = 10®, as 
shown in Fig. ^. The location of the minima have changed 
to ^ = 10^-^ and ^ = lO''. Table | also shows the resuhs of 
these provisional fits. A low reflected fraction seems to be 
preferred for an unionized disc although a fraction closer to 
1 is acceptable in the case of the disc being ionized. The val- 
ues of the photon indices F are consistent with those found 
by Gierliriski et al. (1999) when they include OSSE 7-ray 
observations. 



3.2 Relativistic blurring 




log f 

Figure 3. (High/soft). Results of fitting the ionized disc model, 
without relativistic blurring or additional spectral components, 
to the ASCA data of Cyg X— 1 in the high/soft state. The solid 
curve shows the reduced value for a reflected fraction of 1, the 
dotted line for a reflected fraction of and the dashed line for 
the refiected fraction being a free parameter. When the reflected 
fraction is a free parameter it has been fitted for at each value of 

i- 

3.2.1 Low /hard state 

We fit the ASCA data of the low/hard state between 4-10 
keV with the ionized disc model convolved with the relativis- 
tic blurring model described above. The reflected fraction is 
fixed to 1 and the iron abundance to twice solar. Without 
relativistic blurring we found there to be two local minima in 
around ^ — 10^ and ^ = 10^ **. We first consider an accre- 
tion disc with the inner radius fixed at rin = 6rg, the outer 
radius fixed at 50rg and the inclination a free parameter. 
This improves the quality of fit significantly around ^ = 10^ 
and only slightly around ^ = lO'''* to such an extent that 
the ^ = 10^ solution is now (only just) favoured. If we now 
allow the inner disc radius to be a free parameter we get a 
better fit for very low values of the ionization parameter but 
almost no change at the best fitting values. These results are 
shown in Fig. M. 

3.2.2 High/soft state 

In contrast to the low/hard state Fig. ^ shows that relativis- 
tic blurring does not significantly improve the quality of the 
best fits, and it is possible to achieve an equally good fit with 
or without blurring, with ionization parameters of ^ = 10^'^ 
and ^ — 10^. In the case of cold reflection, however, blurring 
does significantly improve the quality of fit. 



4 FITS TO THE ASCA, CINGA AND EXOSAT 
DATA 

When we fit the ASCA GIS data we find that a satisfactory 
fit cannot be achieved with a disc that has a single ion- 
ization parameter. This is not too surprising since the the 
models is extremely simplified, and the GIS data have a bet- 
ter spectral resolution that the SIS data at higher energies, 
and provide a more stringent test for theoretical models. A 
model in which the disc is divided into two regions, each 
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Figure 4. (Low/hard). Results of fitting the relativistically 
blurred ionized disc model to the 93 Nov ASCA data of Cyg X— 1 
in the low/hard state. The refiected fraction has been fixed at 1. 
The solid curve shows the fit without relativistic blurring, and 
the dashed line shows the improvement with blurring from an 
accretion disc spanning 6 — 50rg. The inclination angle is a free 
parameter, fitted for at each value of 



ASCA CIS 
-Low / hard 



ASCA SIS 
41igh / soft 



Energy (keV) 



10 



Figure 6. Typical ratio plots of good fits to the ASCA data of 
Cyg X— 1 in both the low/hard and high/soft state. 




Figure 5. (High/Soft). Results of fitting the relativistically 
blurred ionized disc model to the 93 Nov ASCA data of Cyg X— 1 
in the high/soft state. The refiected fraction has been fixed at 1. 
The solid curve shows the fit without relativistic blurring, and the 
dashed curve shows the change when blurring from an accretion 
disc spanning 6 — 50rg is included. 



with a diflterent ionization parameter, does provide a satis- 
factory fit to the data. This is not a perfect model, but it 
is more realistic than one in which there is a single ioniza- 
tion parameter. The details of the model are discussed in 
section 2.5. To fit the data with such a model we initially 
allow rtrans aud the relative normalizations of the inner and 



outer components to be free parameters. Fitting this model 
to the data gives the approximate relative normalization of 
the inner and outer components allowing the corresponding 
^trans to be computcd. rtranss IS thcu frozcu at this value and 
the data are fit again. We check the relative fluxes of the two 
components are still consistent with the rtrans. This ensures 
that each of the components is blurred appropriately. 

The results of fitting this model to the ASCA SIS and 
GIS, Ginga and EXOSAT data are summarized in table ^ 
The reduced values of the fits to the ASCA GIS data 
are comparable to those found by Ebisawa et al. (1996) and 
Done & Zycki (1999). Overall this model is able to provide 





High / soft 










Low / hard 


flnn 





Energy (keV) 

Figure 7. Two models that provide good fits to the ASCA data 
of Cyg X-1 in the high/soft and low/hard states. 



a good fit to the available data. The preferred inclination 
angles tend to be small which is in reasonable agreement 
with the optical studies of Gies & Bolton (1986) in which 
the most probable inclination angle is 28° — 38°. 

The low/hard state may be characterized by an ac- 
cretion disk in which the inner region within ~ 13 
Schwarzschild radii is highly ionized, f ~ 10*, and outside 
of which the disk is less highly ionized, ^ ~ 10^. In the 
high/soft state we have only fitted one data set, but that is 
consistent with a relatively cool accretion disk. In all cases 
the accretion disk extends in to the innermost stable orbit 
and has a reflection fraction of unity. Fig. ^ shows typical 
fits to the ASCA GIS data of the low/hard state and ASCA 
SIS data of the high/soft state. The corresponding models 
are shown in Fig. ^ 

Note that in these fits we have not added any contribu- 
tion to the spectrum from thermal disk emission (which will 
have a small effect on these fits since we are fitting above 4 
keV), or any contribution to the iron line from reflection off 
the companion star. This would consist of a narrow iron line 
at 6.4 keV of low equivalent width. Either of these measures 
could improve the quality of fit listed in table ^ 
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/ 
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13 
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57 
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80 
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1.50 
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20 


81 


/ 


67 
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1.44 


4 3+0-2 
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1 7+0.2 
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53 
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4 7+0.3 
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21 


1 7+0.2 
^■'-0.2 
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/ 
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E14 


1.61 


4 3+0-2 
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20 


64 


/ 


75 


b 


2.46 
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6 


9 n+0.2 
^■^-0.2 


13 


36 


/ 


36 



° power law photon index F 
^ 5.0 is the upper limit for this parameter 
accretion disc inclination angle 



Table 3. Fit statistics for our model of a complete (6rg — > lOOOrj) relativistically smeared ionized accretion disc with a reflected fraction 
of unity. The low/hard state data are above the separating line, the high/soft state data below. 



5 TRANSITION REGION 



So far we have assumed a sharp transition between the high 
ionization inner zone and the lower ionization outer zone. 
It is important to try and quantify the sharpness of this 
transition. Some models such as those of Ross & Fabian 
(1993) predict that as the ionization parameter decreases 
outwards there should be a zone in which H-like and He-like 
iron fluorescence is observed. This has twice the fluorescence 
yield of neutral iron, and produces strong lines at 6.67 keV 
and 6.97 keV. Other models, such at those of Nayakshin et 
al. (2000) and Nayakshin (2000) predict the strength of the 
iron line to decrease monotonically with radius, i.e. it is not 
necessary to have an annulus in which strong H-like and 
He-like iron fluorescence is produced. 

To investigate whether an intermediate zone in which H- 
like and He-like iron are produced is consistent with the data 
we modified the ionized disc model as follows. An extra zone 
was introduced of width Ar, with an ionization parameter 
fixed so that the refiected spectrum had H-like and He-like 
iron emission, and the normalization fixed so that the ratio 
of model fiux in 6rg rtrans to that in rtrans nrans + Ar 
is equal to that expected for a standard accretion disc. Fits 
were then performed to see how the best fitting values 
changed with Ar. For data set 3 {ASCA GIS) the results 
are shown in Fig. ^ 

The introduction of a zone of H-like and He-like iron 
makes the fit progressively worse as the width of the zone 
increases. The overall best flt is achieved without such a 
zone, i.e. with a sharp transition. With the present quality 
of data, however, acceptable fits are obtained with a zone 
of reasonable width, ~ few rg, and we cannot rule out a 
gradual change in ionization parameter. 



Width of "hot" Fe zone (r^) 

Figure 8. Change in of flt to data set 3 {ASCA GIS) with 
the introduction of a transition zone of width Ar in which the 
reflected spectrum has H-like and He-like iron lines. The best 
overall fit is achieved with a sharp transition (i.e. Ar = 0) but the 
data are not of sufficient quality to rule out a gradual transition. 



6 DISCUSSION 

We have shown that a relativistically blurred accretion disc 
with two ionization states is able to flt the X-ray data of 
Cyg X-1. The disc may extend in to the innermost stable 
orbit around the black hole, and have a reflected fraction of 
1. We present this model as an alternative to those in which, 
in the low/hard state, the outer accretion disc is not ionized 
and extends in to a few 10s of Schwarzschild radii, within 
which the accretion flow becomes geometrically thick, opti- 
cally thin and radiatively inefficient. One attraction of our 
model is that the accretion flow need not change qualita- 
tively between the two spectral states of Cyg X-1; it just 
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state "5 ''r 

Low/hard High > 0.4 Low Low 
High/soft Low < 0.1 High High 

" ionization parameter 

^ fraction of accretion power dissipated in the corona 
inferred reflected fraction 
observed photon index 

Table 4. Model predictions; see discussion 

requires a change in the surface density of the disc. Fig. ^ 
shows two representative fits to the data. 

It is also clear that, in the low/hard state, relativistic 
blurring of the disc spectra results in a significant improve- 
ment in the quality of fit. 

It should be noted that our model is simplified in that 
the accretion disc is assumed to have two ionization zones. 
With the present data we are unable to investigate the radial 
(or indeed azimuthal) dependence of the ionization param- 
eter in sufficient detail. 

6.1 State changes in Cyg X— 1 

We now discuss a possible accretion disc model to explain 
the differing spectral properties of Cyg X-1 in the low/hard 
and high/soft states. Most of the time Cyg X-1 is in the 
low/hard state and it makes occasional transitions to the 
high/soft state. The mechanism for this transition is at 
present poorly understood. We propose a model in which 
the state change is due to a change in the ratio, /c, of ac- 
cretion power dissipated in the corona to that dissipated in 
the disc. 

If most of the accretion power is dissipated in the disc, 
/c ;$ 0.1, we expect significant thermal blackbody emission 
from the disc. There are then many UV seed photons to be 
inverse Compton scattered in the corona and the X-ray spec- 
trum will be relatively soft. If on the other hand significant 
power is dissipated in the corona, /c > 0.4, we expect the 
blackbody emission to be reduced. The seed-photon-starved 
corona will produce a harder spectrum. A combination of 
factors may then give rise to the observed behaviour of the 
ionization parameter. Firstly, a larger fraction of accretion 
power being dissipated in the corona provides increased hard 
X-ray illumination of the disc, and hence a larger value of ^ 
(note that the bolometric luminosity remains constant with 
changing /c but the ratio of hard to soft flux does change). 
Secondly, the inverse Compton temperature is much higher 
for a hard illuminating spectrum than a soft one, and if 
the outermost layers of the disc heat up to the Compton 
temperature (e.g. see Rozahska 1999; Nayakshin, Kazanas 
& Kallman 2000) this will be important in determining the 
oberved value of ^. A final possibility is that the density 
of the surface layers of the disc changes between spectral 
states, with the low/hard state having a lower surface den- 
sity (giving rise to a higher value of ^) than the high/soft 
state. Table ^ summarizes the expected properties of this 
simple model. 

We also note that the predictions of changing spectral 
slope are consistent with the calculations of Coppi (1999) in 
which an increase in the soft fiux coupled with a reduction 



of the hard flux leads to an increase in the photon index of 
the continuum, i.e. a softer spectrum. 

Our simple model is also consistent with the observa- 
tional results of Zdziarski, Lubinski & Smith (1999) whereby 
the inferred (cold) reflected fraction is correlated with the 
observed photon index. The correlation is almost linear 
ranging from a reflected fraction of with a photon index of 
r = 1.6 up to a reflected fraction of 2 with a photon index 
of r = 2.2, with some scatter. This correlation is also seen 
within individual objects whose photon indices and reflected 
fractions change with time. We predict that a softer spec- 
trum (high photon index) should have a lower value of the 
ionization parameter and hence most of the disc may pro- 
duce a cold line resulting in a large inferred reflected frac- 
tion. In the case of a harder spectrum (low photon index) we 
predict a higher value of the ionization parameter and the 
spectrum will be that of an ionized disc as discussed above. 
This may result in the inferred reflected fraction being lower. 
Another model that is able to reproduce this correlation is 
that of Beloborodov (1999) in which the X-ray sources move 
away from the disc at mildly relativistic speeds. 

Further detailed calculations are, of course, required to 
investigate this model properly. It is not clear how Cyg X- 
1, and other BHC, remain 'locked' into their spectral states 
when the dynamical timescale is so short. 

6.2 Future observations and theory 

Future observations with higher spectral resolution will pro- 
vide much more stringent tests of the accretion disc mod- 
els for Cyg X-1 and other BHC. At present, even with the 
ASCA data, we are unable to place strong constraints on 
the parameter space of our model (e.g. inner disc radius, 
reflected fraction) as our statistics are not good enough. 

This model is also applicable to other black hole candi- 
dates and we hope to investigate these in the future. More 
complete models that account for the detailed ionization 
structure of the accretion disc as a function of radius need 
to be computed. 
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APPENDIX A: BLURRING MODEL SPECTRA 

The disc spectra are convolved with the kernel of the 
DISKLINE model (Fabian et al. 1989). In order to do this 
our disc spectrum must be evaluated outside the range over 
which we are fitting the data. This is because, for example, 
higher energy photons may be redshifted into the band we 
are fitting. If the response matrix of the detector is such 
that the set of channels that correspond to a higher energy 
photon is distinct from the set of channels channels that cor- 
respond to a photon in the energy range being fitted then 
the model will not be evaluated at that higher energy and 
hence that contribution will be lost. This leads to the model 
spectrum dropping off at the edges of the range being fitted. 

A solution to this is to firstly combine the response and 
auxiliary response matrices, and to replace the auxiliary re- 
sponse matrix with one in which the effective area of each 
energy bin is unity. The extend command may then be 
used to enlarge the energy range over which the model is 
evaluated. This solves the problem. 



